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Abstract 


A near-optimal guidance law for the ascent trajectory from earth surface to earth orbit of a 
fully reusable single-stage-to-orbit pure rocket launch vehicle is derived. Of interest are both the 
optimal operation of the propulsion system and the optimal flight path. A methodology is 
developed to investigate the optimal throttle switching of dual-fuel engines. The method is based 
on selecting propulsion system modes and parameters that maximize a certain performance 
function. This function is derived from consideration of the energy-state model of the aircraft 
equations of motion. Because the density of liquid hydrogen is relatively low, the sensitivity of 
perturbations in volume need to be taken into consideration as well as weight sensitivity. The 
cost functional is a weighted sum of fuel mass and volume; the weighting factor is chosen to 
mini mise vehicle empty weight for a given payload mass and volume in orbit. 



Nomenclature 

lh 2 

LOX 

= liquid hydrogen 
= liquid oxygen 

D 

= drag, lb 

M 

= Mach number 

E 

= total mechanical energy per unit 

Mo- 

= transition Mach number 


weight, ft 

Mr* 

= optimal transition Mach number 

8 

= gravitational acceleration on the 

R 

= radius of the earth, ft 

earth surface, ft/sec 2 

T 

= thrust, lb 

h 

= altitude, ft 

Tv 

= magnitude of thrust component 

Isp 

= specific impulse, sec 


along velocity vector, lb 

K 

= weighting parameter, lb/ft 3 

t 

= time, sec 

k r 

= value of K for minimum empty 

V 

= speed, fps 


weight, lb/ft 3 

v P 

= propellant volume, ft 3 


L = lift, lb 
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w 

= aircraft earth surface weight 
(mass), lb 

Wp 

= propellant earth surface weight, lb 

<f> 

= cost functional 

t]hr 

= mass flow ratio of liquid 
hydrogen to liquid hydrocarbon 

T]oh 

= mass flow ratio of liquid oxygen 
to liquid hydrogen 

f]OR 

= mass flow ratio of liquid oxygen 
to liquid hydrocarbon 

P 

= net propellant density, lb/ft 3 


Subscripts 

DF = dual-fuel mode 

E = empty 

f - final value 

H - liquid hydrogen 

LO = lift-off 

O = liquid oxygen 

R = liquid hydrocarbon 

SF = single-fuel mode 

0 = initial value 


Introduction 

Current studies of single-stage-to-orbit 
(SSTO) launch vehicles are focused on all- 
rocket propulsion systems 1,2 (Fig. 1). One 
feature of a SSTO vehicle is its low payload- 
to-gross weight fraction. This means that 
vehicle performance is extremely sensitive to 
perturbations in vehicle design and 
operation. In particular, it is essential to 
“optimize” the flight path and the operation 
of the propulsion system to the extent 
possible in order to attain adequate mission 
performance, and to do this for every 
competing design under consideration. 



Figure 1. Illustration of the SSTO rocket servicing 
the Hubble Space Telescope. 


Dual-fuel (tripropellent) systems that 
initially bum both kerosene and hydrogen as 
fuel and later switch to soley LH 2 , might 
enable attainment of high mass fractions, 
principally because of their greater average 
propellant density and the greater thrust-to- 
weight ratio of LOX-kerosene engines. 3 * 6 
The advantages of hydrocarbon fuel are 
important early in the ascent trajectory, 
where vehicle weight is high, and its use may 
be expected to lead to reduced vehicle size 
and weight. Because LH 2 is also typically 
needed for cooling purposes, in the early 
portion of the trajectory both fuels usually 
must be burned simultaneously. Later in the 
ascent, when vehicle weight is lower, specific 
impulse is the key parameter, indicating 
single-fuel LH 2 use. 

Two recent papers 5,6 have considered 
the optimization of dual-fuel SSTO vehicles. 
Included in the studies was a determination 
of Mj, the Mach Number at which to 
transition from dual-fuel mode to LH 2 
operation in order to minimize vehicle empty 
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weight. Both of these references treat M a as 
an external design variable. 

In this paper, a guidance algorithm is 
developed that determines whether dual-fuel 
or single-fuel operation is superior as an 
integral part of the trajectory integration. 

This approach saves a substantial number of 
iterations of a computer design code by 
reducing the number of design variables, and 
hence the number of design iterations 
required in a vehicle optimization study. 
Further, the guidance law will be directly 
useable as part of a real-time, on-board, 
propulsion control system. 

The basis of the guidance law is the 
energy-state dynamic model. The key idea is 
to introduce the total mechanical energy as a 
state variable, and then to neglect all other 
dynamics. When flight path optimization is 
done with this model, simple rules for the 
optimal path and for the optimal operation of 
the propulsion system are obtained. This 
dynamic model has been used successfully 
many times to obtain effective guidance laws 
for a wide variety of aircraft and missions 
(see Ref. 7 and the references therein for a 
review of this work). The energy-state 
approach is particularly suitable for launch 
vehicles because efficient energy 
accumulation (or equivalently maximizing 
“total A V”) is the primary trajectory 
optimization goal. 

In a series of papers 7 ' 9 we have used 
energy-state methods to develop algorithms 
for ascent trajectory optimization and 
optimal operation of single-fuel multi-mode 
propulsion systems. In particular, the 
operation of propulsion systems with two 
separate engines, air-breathing and rocket, 
was investigated 9 . The present paper 
extends those methods to the dual-fuel case. 
The main goal is to determine M a and to 
investigate optimal trajectories. 


In the numerical results, vehicle 
performance is computed using the NASA 
Ames hypersonic vehicle synthesis code 
(HAVOC) 10 . HAVOC integrates geometry, 
aerodynamics, propulsion, structures, 
weights, and other computations to produce 
point designs for a wide variety of launch 
vehicles. It is capable of iteratively 
determining “closed” vehicles , that is, 
designs which meet specified payload mass 
and volume requirements for a specified 
mission. Although the trajectory guidance 
law is based on the energy-state model, the 
trajectory integration in HAVOC uses a 
point mass model, including the effects of 
earth rotation, earth curvature, and variable 
gravity. 


Optimization Function 

The energy-state model is obtained by 
using the total mechanical energy per unit 
weight as the state 
variable 7 ' 9 : 

E=P (1) 


• T 

W=-— ( 2 ) 

*SP 


where 


E= 


bR 

R+h 



( 3 ) 


and 


P= ~{Ty- D) (4) 
w 

and where the drag is evaluated at the lift 
required for equilibrium of forces 
perpendicular to the flight path. 


305 



( 9 ) 


For a SSTO mission, what is desired is a 
trajectory that gives the minimum empty 
weight vehicle to put a given payload mass 
and volume in orbit. Because the density of 
liquid hydrogen is low, the sensitivity of 
perturbations in volume need to be taken 
into consideration as well as mass sensitivity, 
and it is therefore necessary to minimize a 
weighted sum of propellant weight and 
volume. Thus we introduce the cost 
functional 


t=W p+ Kv f (5) 

where K e [O,®) is a weighting parameter to 
be chosen later. 

The quantity to be minimized for a given 
energy gain is 

S= (</* = '( }dt= f-dE (6) 

where Eq. (1) was used. It is assumed that 

^ >0, P> 0 and that E is monotonically 
increasing. If the propellant density is 
P = W P I v r , then from Eqs. (2) and (5), and 

• m 

using W P = -W, 


• • W. T 

<t> = W p +K- L = 




' ^ 
1+— 
P) 


( 7 ) 


For convenience, we choose to invert the 
integrand in Eq. (6) and maximize; from Eqs. 
(1), (6), and (7), the quantity to be 
maximized is 



( 8 ) 


where 


W 7 1+— 

V PJ 

The guidance algorithm then consists of 
selecting propulsion system and trajectory 
parameters that maximize the function F as 
given by Eq. (9) at each energy level along 
the trajectory, subject to any relevant 
constraints. 

For vehicles capable of either dual- or 
single-fuel operation, the densities to be used 
in Eq. (9) are 


Pdf = 


PrPoPjX 1 + *1 on 
( PqPh + 7 orPrPh + *1 hrPrPo) 


( 10 ) 


PSF = 


PoPjM^ 7 qj /) 

{Po + VohPm) 


( 11 ) 


SSTO vehicles are typically subject to 
dynamic pressure constraints and a maximum 
tangential acceleration limit. This latter limit, 
nominally 3 times the earth surface 
gravitational acceleration, is met by engine 
throttling. It may happen that the limit 
affects dual-fuel operation but not single-fuel 
operation at a point along the trajectory. All 
these constraints are accounted for in the 
guidance algorithm. 


In addition to being a useful tool in 
preliminary design studies, the guidance 
algorithm should be ideal for use in an on- 
board real-time control system because : (1) 
it is fully nonlinear and models all of the 
vehicle's significant nonlinearities, (2) it is 
algebraic and thus does not rely on 
potentially unstable numerical integrations, 
and (3) it depends directly on easily 
measured vehicle states and parameters. 
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Numerical Results 

All numerical examples will be based on 
an SSTO rocket with a delta winged-body 
configuration 2 (Fig. 1). The three 
propellants (hydrocarbon fuel, LH 2 , and 
liquid oxygen) are stored in three separate 
internal tanks. The vehicle takes off 
vertically and lands horizontally. The first 
results to be presented use a fixed trajectory 
commonly used for SSTO rockets. 

As a first step, the best transition Mach 
number, M/, is determined by treating this 
parameter as a single external design 
variable, as was done in Refs. 5 and 6. The 
results are shown in Fig. 2, which plots gross 
lift-off weight ( W L o) and empty weight ( W E ) 
as a function of Mu It is seen that both 
minimum Wlo and We are obtained at about 
M* = 9.0, and that the weight savings at 
Mu are substantial relative to low values of 
Mu All of the data points on Fig. 2 are for 
closed vehicles and hence several design 
iterations are necessary for each point. 



Figure 2. Effect of transition Mach number on 
vehicle gross take-off weight and empty weight. 

Before applying the developed guidance 
law to this problem, the best value of K must 
be determined. This is done by computing 


closed vehicles for a range of values of K 
(Fig. 3). It is evident that a value of K = 4 
lb/ft 3 , denoted hereafter by K\ gives very 
nearly a minimum of both empty weight and 
gross lift-off weight, and this value will be 
used throughout the rest of the paper This 
value of K represents a factor of over 1 0 in 
weighting the cost functional in favor of 
propellant mass (a value of K= p would 
signify equal weighting of propellant mass 
and volume.) The figure shows that the use 
of the optimally weighted cost functional 
saves 1.7% in empty weight and 1% in gross 
lift-off weight, relative to minimizing 
propellant weight only. 



It is of interest to compare these results 
with the equivalent results for an air- 
breathing launch vehicle, as shown in Fig 2 
of Ref. 9. For the airbreather, the best value 
of A: is also around 4, but the empty weight 
reduction relative to minimizing propellant 
weight only is much larger, at 4.9%; this is of 
course because all of the airbreather 
propellant is low-density LH 2 , and therefore 
this vehicle is more sensitive to volume 
perturbations. 

Figure 4 plots the function F along the 
fixed trajectory. Whichever mode of 
operation, dual-fuel or single-fuel, that gives 
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the highest value of F at a given speed 
should be the one selected at that speed. 

The figure shows that from lift-off to M= 
9.0, the dual-fuel mode is superior, and 
above this speed the single-fuel mode is best. 
This value of M- = 9.0 agrees with the value 
determined by treating M a as a design 
variable. Fig. 2, thus validating the guidance 
law. The value of M-* as determined in Ref. 

5 was in the range 8.6 - 8.9, and for Ref. 6 it 
was in the range 7.3 - 7.4. 


rs Mach No. @K= 4 



Figure 4. Cost function histories for both propulsion 
modes. 

The relative distance between the two 
curves on Fig. 4 provides an assessment of 
the difference in performance between the 
two modes at a given Mach number. It is 
seen that both modes give substantially the 
same performance between M= 7 and M= 
11. This relative insensitivity to 
characteristic of a design variable near its 
optimal value was also observed in Ref. 5. 
The use of single fuel LH 2 mode becomes 
increasingly advantageous as Mach number 
increases past 1 1 . 

Fig. 5 shows the effect of K on M- As 
expected, the higher the value of K, the 
higher the premium on minimizing fuel 
volume, and thus the larger the transition 
Mach number, M 



Figure 5. Effect on K on Mtr. 


The function F was also used to 
optimize the ascent trajectory (Fig. 6). As 
compared with the fixed trajectory, the near- 
optimal one has increased dynamic pressure, 
especially in the initial dual-fuel mode. 



Fig. 7 shows the cost function surface 
within the flight envelope. The optimal 
energy climb path begins at zero altitude and 
then encounters the maximum dynamic 
pressure boundary and begins to climb along 
this boundary up to Mach 3.2. After this 
point the trajectory increases in altitude and 
is not constrained by either the upper or 
lower dynamic pressure limits. At Mach 9.6, 
the tripropellent system switches to single- 
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fuel mode, and the altitude transits to a 
higher level. 



Figure 7. Cost function surface and optimzed flight 
path. 


The transition Mach number along the 
optimal path is equal to 9.6, as can be seen 
from Fig. 8. 



Figure 8. Cost function histories of optimal energy 
climb path. 

Fig. 9 shows the weight ratios at each 
energy level. It is seen that at orbit energy 
level the optimal flight path gives a small but 
significant savings of fuel relative to the fixed 
path. The near-optimal trajectory consumed 
less fuel in the amount of 0.9% of W LO than 


did the fixed, almost all the difference 
occurring in dual-fuel mode. 



Figure 9. Weight ratios at energy levels. 


Concluding Remarks 

A simple guidance law for operation of 
dual-fuel SSTO launch vehicles has been 
developed and used to determine the optimal 
value of the transition Mach Number from 
dual-fuel to single-fuel. For the example 
considered, the optimal transition Mach 
Number was 9.0 along a fixed trajectory. 
Along an optimal trajectory, the best 
transition Mach number was 9.6; the optimal 
trajectory had higher dynamic pressure than 
the fixed, particularly in dual-fuel mode. 

In the future, the guidance method 
described in this paper easily could be 
extended to optimize other propulsion 
system parameters, such as flow rates of 
individual propellants in multi-propellant 
engines. Because the guidance algorithm is 
internal to the trajectory optimization 
routine, its use will save many iterations of a 
preliminary design computer code relative to 
treating these parameters as external design 
variables. The guidance law is highly 
accurate, robust, and simple to implement, 
also making it ideal for use in a real-time on- 
board control system for a SSTO launch 
vehicles. 
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Introduction 

Application of optimal control theory in the form of the maximum principle to aircraft 
trajectory optimization problems generally results in a two-point (2PBVP). The order of 
this problem is double the number of state variables and the equations are always “half 
unstable.” Many schemes have been developed to numerically solve this difficult class of 
problem, but all are unsuitable in a vehicle synthesis code. Not only are they 
computationally expensive, but they are non-robust and not user-friendly. 

What is needed in a vehicle synthesis code is a method that optimizes the trajectory in 
one pass, that is as an integral part of the trajectory integration. The method must also be 
robust and it should be easy to use and to interpret physically. The key to achieving this is 
to use judicious approximations to reduce the functional optimization problem to a 
function one. 

In this report, two approximation techniques are reviewed and developed. The first is 
the use of the Energy State Approximation (ESA). This well-known technique substitutes 
the total mechanical energy for the speed as a state variable, and then neglects the altitude 
and flight path dynamics relative to the energy dynamics. The second technique is the use 
of Singular Perturbation Theory (SPT) to time-scale decouple equations of motion. These 
two techniques are related, and in fact, the ESA may be viewed as an example of SPT 
methods. 


1 



Trajectory Optimization: the Maximum Principle 


The equations of motion of aircraft flight, no matter what the assumptions (see 
Appendix A), are of state equation form: 

x = f(x,u) 


where X_ s 91" is the state and U_ eUCSK” is the control vector. Suitable boundary 
conditions on the state vector components are prescribed. It is desired to find the 
components of U along the trajectory so that 

J - £' fo(X, U)dt 

is minimized. It is assumed that time is free. The necessary conditions for optimal control 
are provided by the Maximum Principle (MP). 


Theorem (the Maximum Principle): Introduce the variational Hamiltonian function 

h = -A+j-t,/;. 

i=\ 

where the components of the adjoint vector, A - , satisfy the differential equations 



d_ H 
dX . 


, / =\, n 


Then, if U is an optimal control, 


(a) Lf - arg max H 
Us U 


(b) 11 = 0 

(c) Transversality conditions (“natural” boundary conditions on the A,) hold. 
Thus, we must solve a 2n dimension 2PBVP in the states and adjoints; exactly n 
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boundary conditions are provided at t = 0 and the other half at t = tf (due to the 
travsversality conditions). The equations are unstable in the sense that if they are 
linearized about a nominal trajectory, one-half of the eigenvalues will have positive real 
parts and the other half negative (unless they are zero). 

A pproximation Techniques 

Methods of reducing the 2PBVP to a simpler problem will now be developed. These 
methods focus on order reduction and are motivated by two simple observations. 

First, we note that if all componets of f , except possibly /}, and the function f Q are 
independent of a specific state variable, say X ; , and the final value of X, is not specified, 
then the corresponding adjoint is always identically zero and the state equation X t = f t 

drops out of the problem (decouples from the other states). To see this, consider the 1 
adjoint equation and its travsversality condition: 



±h-i, 

d X. ' 


A i (tf ) “ 0 


The only solution to this problem for any finite value of d f i fdX i , is A, = 0. 


Second, we note that if there is only one state equation, then the necessary conditions 
can be used to eliminate the adjoint variable and thus the problem reduces from a 
functional optimization problem to a function one. To see this, consider 


X = f(X,U) 


J = £' fo(X,U)dt 


We have 


H = -/o + A / 
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■ x .**-xlL 

d X d X 

Applying the MP (assuming for the moment unbounded optimal control exists) 


Yi = -fo + X f = 0 

d H _ d fo d f _ 
d U ~ d U +A d U 

Eliminating X gives 

d fo . df r A 
f + fo - 0 

d U d U 

for the optimal control. Alternatively, a direct approach may be used: 

Thus ( fo/f ) is to be minimized with respect to U at constant X, this leads directly to the 
equation for optimal control derived just above from the MP. 


SPT provides an organized, mathematical way to view order reduction of differential 
equations. Consider the initial value system 

X = f(X, Y) X(0)=X 0 

eY = g(X,Y) (0 )=Y 0 

where e is a “small” parameter. 

Since e is small, an approximate system may be expected to be 
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X,=f(X„Y r ) 
0 = g(X„Y,) 


It can be proved that under certain conditions, the solution of this problem is a good 
approximation to the solution of the original problem, except near t — 0 because the 
boundary condition K(o) = Y 0 will be generally violated these. 

The problem is that Y undergoes a rapid transition from its boundary condition to the 
approximate solution at t = 0 . To analyze this motion, the time scale is stretched by 
T = t / e. The resulting equations are called the boundary layer equations 


dX 

dT 


ef(X,Y ) 


dY 

dT 


S(X,Y) 


Setting e = 0 to approximate these equations 


— = 0 => X = const = X. 
dT ' 


dT 


= g{X>,Y b ) 


The solution to this equation approximates the desired solution near t = 0 . There are 
matching techniques to combine these two solutions to give an over-all approximation, if 
desired. The key observation is that a second order system has been replaced by two first 
order systems, and each of them reduces to a function optimization problem. 
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The SPT provides a convenient way to look at the energy state approximation. Define 
the aircraft energy per unit weight by 

E = h + — V 2 

2g 

Differentiate and use the state equations in the Appendix 

t. 

g Mg 

Where T v is the component of thrust along V_ and P is the specific excess power. Note 
that this equation is valid for all three sets of equations given in the Appendix. 

Now replace V by E as state variable and use the observation that h and y are 
capable of rapid change relative to E . This motivates writing 


E = P 


eh =■•• 
e y =••• 

Setting e = 0 then gives an order reduction of the equations of motion by two. We will 
use this approximation, the ESA, throughout. 

This approximation has a long history of successful application in a wide variety of flight 
trajectory problems. The main drawback is that the variables h and y may now jump 
instantaneously at points along the trajectory, as well as at the boundaries. These jumps 
could be accounted for by boundary layer analysis, but this is not done in this report. 
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ESA equations of motion for the three cases of interest are given in the Appendix, 
these equations will now be used as the basis for discussing specific trajectory 
optimization problems. 

Minimum Time/Fuel to Climb 

Starting from equations (1)' 


m = -CT 


m(0) = m 0 


X = V 

E = P E(0) = E 0 ,E(t f ) = E f 

L = mg 

J = \{K X + K 2 CT)dt 

Here, the system functions and boundary conditions do not depend on X and thus the 
equation X - V drops out of the problem: 

m = -CT m(0) = m 0 

E = P E(0) = E 0 ,E(t f ) = E f 


/ = J (AT, +K 2 CT)dt 

with P evaluated at L = mg . The system functions and boundary conditions now 
depend on both Eand M and hence neither state equation uncouples. Thus the MP must 
be applied and a 2PBVP solved. 
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To reduce the problem to one of function optimization, SPT is used to further reduce 
the system. 

m = -CT 
eE = P 

Setting e = 0 gives a single state equation 

m = -CT m(0) = m 0 

T = D, L = W 

The optimization problem is now (h and E are controls) 

/ = (/:,+ K 2 CT) t f 

With the obvious trivial solution t f - 0 . Also, because the system functions do not 
depend on m and m(f / ) is free, X u =0 . 

The boundary layer system for this problem with e - 0 is simply 

E = P 

with m = const, so that 

/ = j(K,+K,cr)^- 

and the solution reduces to 
max/' p 

h + KjCT ) E=coasJ 
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assuming that + ^ * s positive and £"is monotonic. This is the well-known 

energy climb path. 

From now on we will assume “slowly varying” mass, that is that m is on a slower time 
scale than E and thus its state equation may be ignored. It is also assumed that the throttle 
is fixed. 

Fixed Ranee 

This problem is the same except that range is fixed 

X = V X{0 )=X 0 ,X(t f ) = X, 


E = P 


E{0) = E 0 ,E(t f ) = E f 


J = \{K ] + K 2 CT) dt 

Here, X x - const 4- 0 so that the X - V state equation does not uncouple, and we have 
a 2PBVP. To effect system order reduction, SPT is used. 

X = V 

e E = P 


H - -AT, - K 2 CT + X x V + X e P 


= 0 

* ax 


X x - const 
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e X 


E ~ 


dH „ d{CT) , d V , <? P 

XV ^ A v - ^ ^ 

d E d E d E E d E 


Setting e = 0 a problem with a single state is obtained 


X = V 


T=D 


H = -AT, -K 2 CT + X x V 


L=W 


Applying the MP: 

Min (K x +K l CT'\ _ K x + K 2 C c T c 

h,E l V J t w " F c 

AT, +K 2 C c T c 

~ y 

v c 

This defines a cruise point, characterized by C c , T c , and V c , in the flight envelope. By 
proper selection of Ki and K 2 , this point can be made to closely approximate minimum 
direct-operating-cost cruise. 

For minimum time (K x =1, K 2 = 0) , the optimum cruise point is given by 
Max 

h,E 

For minimum fuel consumption (A - , = 0, K 2 = l) it is 



Max f V Max 

(V{Vof 1 

h, EKCT J ~ h, E 

l c J 


which is the classic Brequet cruise point. 


10 



The boundary layer with e = 0 is 


E = P 


H = -K^ - K 2 CT + X x V + Z E ? 
so that the optimal climb flight path is given by 


Max/'’ P ^ 

h [k, +K 2 CT-X x Vj 


E= const 


with X x as given above. 


Maximum Turning With No Thrust 


We start with (2)' with T=0. 


X - V cos x 


Y = V sin x 


E = P = - 


VD 

mg 


X = 


L sin <f> 
mV 


L cos <f> = mg 
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with 


J = " J X dt = -J 


L sin ^ 
mV 


In this case the X, Kand £ equations all uncouple. Changing to more convenient 
variables: 


where 


E = -V(b + Cco 2 ) 


J - -^(o fdt 


(O - 


tan <j> 
tan <f> M 


- 1 < <y < +1 


<f> M = sec l [min(C i limit, load factor limit)] 


f = 


g tan <f> 


M 


Thus 


'■! 


cof 


v[b + Cco 2 ) 


dE 


and the optimal controls are given by 


Max 

h,o> 


! (Of 


v(b+cw 2 ) 


I E=coost 
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which leads to 


d_ 

dco 


( )-0=>.-±g 



(assuming right hand turn) 


dh 


() = 0 : 


V d h V) 


( B + Ca 2 )-/vf— +—c o 2 
' ' \d h dh 


= 0 


where 

g Do±DlO C Dj^lSl 
Mg ' Mg’ 


so that 




The search for the optimum h probably should be done numerically. 

Next, consider the same problem but using (3)' . Now, the X equation uncouples 
but the ' additions do not. The coupling is of two types. First, through the 

Coriolis terms, which are relatively small and can be ignored. Second, through the 
centripetal terms, which are large at the start of descent trajectories from orbit. 
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There are two ways to deal with this problem. First, the Coriolis and centripetal 
terms are ignored. This is justified because what is really sought is turning ability due to 
banking and these terms mask this. Second, the E and x terms may be decoupled using 
SPT. 

Maximum Cross-Ranee 

Next consider, using (2)' 


/ = -JY<# =-jVsinxdt 


X - V cos^f 

= Vsinx 

E = P 

Ls\nd> 

X = v 

mV 

Lcostf) = mg 

As before, the Xand Y equations uncouple but now the % equation does not. To 
reduce this to a function optimization problem, further time-scale separation is required. 
Putting x on a slower time scale than E gives the solution V = 0 and k x - 0 . 

t 

Using equations (3) results in the same problems as for maximum turning. 
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Appendix 

Equations of Motion 

The following are the aircraft point-mass equations of motion under various 
approximations. 

(1) Flight in a vertical plane over a flat, non-rotating earth; no winds aloft and 
thrust aligned with velocity. 


737 = -CT 


X - V cosy 

h = V s\x\y 

■ _ T - D - mgsiny 
m 

. _ L - mg cosy 
r ~ mV 

(2) 3-D flight, otherwise the same as (1). 
m = -CT 
X = V cosy sin x 

= V cosy cos^y 
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h - siny 


• _ T-D- mg s\ny 

m 


L s\n<f> 
mV cosy 


r = 


L cos <j>- mg cosy 
mV 


(3) 3-D flight over a spherical, rotating earth; no winds aloft, thrust not aligned 
with velocity, terms in the square of the earth rotation ignored. 


T 

m = 

& SP 


V = 


T cos P cos (a +£)-£> 

m 


- gsmy 


f T cos /? sin(a + £) + L 




r = 


■ 


mV 


g 

I COS <t>~— cosy 


V Y 

+ — cos/ + 2ru cos /cos— 


'T cos/?sin(ar +£') + Z > 
^ mV cosy J 


. , V Y 

sin 0 cos/ cos / tan — 


_ f . Y . Y] 

+2<u|^tan y sin / cos— - sin — J 
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h = Vsiny 


^ _ VRcosy cos^ 

X ~ Y 

rcos — 

R 

y. _ VRcosy sin^ 

r 

The following are the energy-state approximations of these equations. 

(O' 

m = -CT 

x = v 

P 

mg 

L- mg 

( 2 )' 

m = -CT 
X - V co sx 

= V s\nx 

E = P 

L sin 6 

X — 

mV 

Lcos<f> - mg 
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(3)' (with T = 0 and m - const) 


E = - 


VD 


m. 


_ L . g V „ Y 

O = cos^- — + — + 2co cos ^ cos — 

mV V r R 


L . , V Y . Y 

sin© cosy tan 2<ysin — 

mV r R R 


X = 


V7?cosy 

F 

/•cos — 

R 


VRsin x 


r 
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